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We present higher-derivative gravities that propagate an arbitrary number of gravi- 
tons of different mass on (A)dS backgrounds. These theories have multiple critical 
points, at which the masses degenerate and the graviton energies are non-negative. 
For six derivatives and higher there are critical points with positive energy. 
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1 Introduction 

Two-derivative Einstein gravity in four dimensions is non-renormalizable. It can be made per- 
turbatively renormalizable by adding four-derivative terms to the Lagrangian [l|,|2j. However, 
the addition of the curvature-squared terms spoils unitarity: around a Minkowski background 
they introduce a massive spin-0 and spin-2 mode. These massive modes have norm opposite 
of the massless spin-2 mode, and thus are ghosts. The spin-0 mode can be eliminated by 
tuning the coefficients of the curvature-squared terms, but the massive spin-2 modes cannot. 

Recently, a consistent four-derivative theory of gravity in three dimensions, called 'New 
Massive Gravity' was introduced in [3j. NMG is ghost-free due to the fact that massless 
gravitons have no propagating degrees of freedom in three dimensions, which makes it pos- 
sible to choose the overall sign of the action such that the massive gravitons have positive 
energy. This, however, is not possible in higher dimensions as there both massive and massless 
gravitons propagate. 

One way around this problem is to perturb around an (A)dS background, instead of 
a Minkowski background. The cosmological constant and the coefficients of the curvature 
squared terms can then be tuned such that the massive modes becomes massless [4j. This 
is known as 'critical gravity' [5j. As the massive modes disappear at the critical point, the 
theory is potentially unitary. 

However, at the critical point the massive modes are replaced by so-called log modes [g, L7j, 
. As it turns out, these log modes are ghosts [9], and must the truncated to restore unitarity. 
As their falloff in the radial AdS coordinate is logarithmic (hence the name) , this may be done 
by imposing certain boundary conditions. 

The resulting theory is then unitary, but, unfortunately, also empty. Namely, at the 
critical point the energy of the massless graviton modes vanishes, together with the mass of 



the Schwarzschild black hole. It was recently argued from a CFT perspective 10J that this 
is essentially due to the fact that critical gravity is of rank two (with the rank being half 
the number of maximum derivatives). Instead, gravity theories of odd rank should not suffer 
from this 'zero-energy-problem'. 

The purpose of this paper is to investigate the criticality conditions for higher-rank theories 
of gravity. It is organized as follows. We first give a non-linear Lagrangian for arbitrary rank 
r, that, on (A)dS backgrounds, propagates one massless and r — 1 massive gravitons, but not 
the scalar ghost mode. Next, we show that the quadratic perturbation of this Lagrangian 
and its linear equations of motion can be concisely written in terms of the so-called Schouten 
operator. This reformulation enables us to calculate the global charges (such as black hole 
masses) and graviton energies for arbitrary rank. From the latter we deduce that the theory 
is critical, i.e. all energies are non-negative, whenever sufficiently enough graviton masses are 
degenerate. In general there will be more than one critical point; hence the name polycritical 
gravities. 



2 Non- linear action 

For a gravity theory of rank r (thus containing at most 2r derivatives) , we would like its 
linear equations of motion to bqll 
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n=0 



□ - 2A - ml) hp, = 0, (la) 

v"V = o, (lb) 

gTh^ = o. (lc) 



This a straightforward generalization of the Fierz-Pauli equations of motion for a single mas- 
sive graviton [11;]. Here however we have r graviton modes with a priori different masses 
mo, mi,...,m r _i. Because these equations of motion should follow from some covariant 
non-linear Lagrangian, one of the masses (say mo) will always be zero. This is due to the 
diffeomorphism invariance of the non- linear theory. We will use the index n = 0, 1, . . . , r — 1 
to indicate all gravitons, and the index i = 1, . . . ,r — 1 for only the massive gravitons. 

In four dimensions and higher, a non- linear Lagrangian whose linearized equations of 
motion are those given above, is 
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R-(d-2)(d-l)A + ^C^ pa [ \j>,U'- [ )('"'■''" 




(2) 



Here C denotes the Weyl tensor, and the coefficients a% = aj(r, d, A, m,j) are functions of the 
rank, the dimension, the cosmological constant, and the graviton masses. For r = 2, this 
action was already written down in [12|. We will give explicit values for the coefficients a« for 
r = 3 below. 

Note that we use the canonical sign for the Einstein-Hilbert term in the above action. 
Flipping its sign is equivalent to changing the overall sign of the action, upon redefining A 
and Oj accordingly. Such a change of sign also changes the sign of the energy of the solutions 
(scc lscction 6ft . and, as noted in the introduction, is particularly important in the d = 3, r = 2 
case [3J. There it is customary to leave the sign of the Einstein-Hilbert term arbitrary. Here, 
however, we have fixed the sign, keeping in mind that we can always flip the overall sign of 
the action in order to choose which modes have positive energy and which negative. 

We have two main reasons for using only Weyl tensors in the higher order terms. Both stem 
from the fact that the Weyl tensor vanishes identically on (A)dS spaces. First, this ensures 
the uniqueness of the (A)dS vacuum. Second, for perturbations around such a background 
the higher order terms do not contribute to the trace of the equations of motion. This comes 
about as follows. 

The full non-linear equations of motion that follow from (|2|) are 

E i»> = --^i^ = G ^ + Y, a * K U = o- (3) 



-gSg- 



Here G^ u is the cosmological Einstein tensor (see Appendix A ) , and K % are the contributions 



from the higher-order terms. Suppressing indices on the Weyl tensor, these contributions have 
the generic form 

K i = 2 ^-n l c + c f|H! _ l g ) c. (4) 



1 See | Appendix A] for our conventions on linearization. 



Thus K l consists of a part that is linear in the Weyl tensor, and part that is quadratic. 
For an (A)dS space, both parts are zero. So on these backgrounds just the Einstein-Hilbert 
contribution of Q survives. This uniquely fixes the background curvature to be A. 

For linear perturbations around (A)dS solutions, the part of K l that is quadratic in the 
Weyl tensor vanishes. The linearized higher-order contributions to E% v come then only from 
the first term on the right-hand side of (HJ), which evaluates to 



[Kl 



V'V'EfC^. 



The linear Weyl tensor C is, just like its non-linear variant, traceless. Upon taking the trace 
of the linear equations of motion, it follows that the linear Ricci scalar vanishes on-shell: 



sTE t 



L 



1 



R L = 0. 



(6) 



As in Einstein gravity, this allows us to impose the transverse traceless gauge [13], i.e. equa- 
tions (|lb|) and (flc|) . for the linear graviton fluctuations hau- Hence the scalar mode h, which 
would otherwise be a ghost, does not propagate. 

In the remainder of this section we will show that the linear equations of motion take the 
form (|la|) . and give explicit values of the Lagrange parameters <Zj for the rank r = 3. The 
linearized equations of motion can be written entirely in terms of Gz v and R by using the 
identities 
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R l 



(7a) 
(7b) 



2d-l 

The former follows from commuting covariant derivatives, while the latter is a consequence 
of the Bianchi identities. Furthermore, in the transverse traceless gauge the linear Einstein 
tensor takes on the form 
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p,v 



(□ - 2A) h 
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(8) 



on-shell, the linear equation of motion E^ v is a polynomial in □ that acts 



And as R L 

on h^ u . We may always choose the parameters at such that it factorizes into the form (fTa 

Indeed, for r 



3 the linear equation of motion becomes 
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(9) 



where the parameter r is given by 



r— i 

n[ 



mf - (d - 2)A 



(10) 



and the squared masses by 



m\ 



ml 



dA 
~2~ 
dA 
^2 



a\ 



+ Vb 



2a 2 
a\ — y/b 



2a 2 



ai + (3d - 4)Aa 2 



d-2 
d-3 



a 2 . 



(11a) 

(lib) 
(lie) 



Inverting the above equations for ao and a\ gives finally 

ai = -^-(mj + ml + dA), (12a) 

t a — 6 

ld-2 . ,, 

„ 2 = + — . (12b) 

The factors of <i— 3 in the denominator indicate that the non- linear Lagrangian ([2]) is only valid 
for d > 4. This makes sense, as the Weyl tensor vanishes identically in three dimensions and 
lower. Similar explicit values for the parameters a% of higher-rank theories can be computed 
along the same lines. While these explicit values are needed for the non-linear action, they 
are not for its quadratic perturbation. As we will see in the next section, the latter can be 
written concisely in closed form using the mass parameters rrn instead of the parameters a^. 
Furthermore, the quadratic Lagrangian will also be valid in three dimensions. 

3 Quadratic action 

Before we set out to calculate the conserved charges and energies of our higher-rank theory, 
it is convenient to rewrite linear equations of motion a bit. We start by rearranging the 
quadratic perturbation of the non-linear Lagrangian ([2]). It is given by 

C 2 = - Vg£, - iv^VV hTa.Lr 1 J C M V- (13) 

We have dropped a total derivative in the Einstein part, and expanded one of the linear 
Weyl tensors in terms of the graviton fluctuations hup • There are more contributions to this 
expansion than V ff V"/i w ', but they drop out because of the contraction with the other Weyl 
tensor. 

Like the linear equations of motion, the quadratic Lagrangian can be written entirely in 
terms of (?£, and R by using the identities ([7J). The resulting expression can be simplified 
further to 

£2 = ~G% (jl 25 + mA o V- (14) 

Here r is given as in (|1U|) . and we have introduced the Schouten operator S. It is defined such 
that when it acts on the graviton fluctuations huu> it gives the linear cosmological Schouten 
tensor: 

S V = S L av . (15) 

The cosmological Schouten tensor S^y is in turn defined such that for vanishing A it reduces 
to the normal Schouten tensor, and that it is zero on (A)dS backgrounds, i.e. S^ u = 0. See 



also Appendix A Surprisingly, the quadratic action (|14p is also valid in three dimensions, 
whereas the non- linear action ([2]) was not. For d = 3, r = 2 and A = it coincides with the 
quadratic action given in [14| . 

Before deriving equations of motion from (|14p , we first list some useful properties of the 
Schouten operator and the Einstein operator Q. The latter is defined in a similar fashion as 
the Schouten operator, 

9 o V = G>- (16) 



For arbitrary symmetric tensors An V and Bn U , we have 



B^SoA 
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SoB^+ li-l^^uB - ^V'V'B 
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ld-2 
2d-l 



VuV„g pa GoA 



pa- 



(17a) 
(17b) 
(17c) 
(17d) 



In the first two lines we have dropped a total derivative while integrating by parts. Note 
that for B^y = Q o A^ v the last term of (|17ap vanishes, and the middle term is the same as 
the commutator (|17d|) . Thanks to this subtle interplay between the Schouten and Einstein 
operators, the linear equation of motion reads 



E 



pv 



5C 2 
5h^ 



'-Qo^S + m^ 
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pv 



(18) 



Upon taking the trace of these equations, we should recover ([6]), that is, R = 0. To see how 
this comes about, we need three additional properties of the Schouten and Einstein operators: 



grgoA ia , = -{d-l)grSoA ia ,, 
ld-2 



g^SoA. 



^ 2d - 1 L 

from which it follows that 

-n V - - ■ d-2 



VT^ -UA-{d- l)KA 



g^QoSoA, 



^v- 2 

A short calculation shows that we indeed recover ( 

'r-l 



AfgoAw. 



r v E% = ±sT (U^- (d-2)AJ go v = (i- £) R l 



0. 



(19a) 
(19b) 
(19c) 

(20) 
(21) 



As the linear Ricci scalar vanishes on-shell, we may go to the transverse-traceless gauge (|lb|) . 
(fTc|) . In this gauge the Schouten and Einstein operators become equal (compare equation 

(ED), 

i 



S o h 



pv 



g oh 



pv 



--(□-2A)V- 



The complete linear equations of motion (|18p can then be written as 



(-1)' 

2r 



JJ(Q-2A- 771^)^ = 0, 

sr hpv = o, 



n=0 



(22) 



(23a) 
(23b) 
(23c) 



with 777-0 = 0. 



4 Conformal invariance 

The overall factor - in the quadratic Lagrangian (|14|) comes from demanding that the Ricci 
scalar in the non-linear action ([2]) has the usual normalization. The advantage of this normal- 
ization is that we recover Einstein gravity upon decoupling the massive gravitons by sending 
their the masses to infinity. We will see later in sections [5] and [6] that also the conserved 
charges and graviton energies reduce to their two-derivative 'Einstein' values in this limit. 

However, an obvious drawback of the overall factor - is that it has poles at the mass 
values 

m 2 i= (d- 2)A. (24) 

One easy way to get rid of the poles is to simply replace the overall factor - by some other 
factor — that has the same mass dimension, but no explicit dependence on m, and thus no 
poles. We can then freely let the masses take the values (|24p . with the drawback that we 
do not recover Einstein gravity upon decoupling the massive gravitons. Another possible 
drawback could be that for the mass values (|24|) the trace of the linear equations of motion 
()2ip vanishes identically, and does not eliminate the scalar mode of the graviton. 

Luckily, the latter does not happen. Instead, for the values (|24|) the linear theory develops 
a conformal invariance. To see how this happens, consider the linear conformal transformation 

Suhfo, = g^v u. (25) 



We would like to know the variation of the equations of motion (|18p under this transformation. 
To this end we first compute the variation of a single Schouten operator, 

S u (S o hfj, v ) = S o (g^ u uj) = — (Ag^ u + V M V„) uj. (26) 

Next, we notice the identities 

go(V fl V v co)=So(V fl V v u)=0. (27) 

Thus for the repeated composition of the Schouten operator only the first term on the right- 
hand side of (|26p is important. The variation of the equations of motion (|18p then becomes 

S U E^ = - ( TT m? - (d - 2)A ] (L (G o V) • (28) 




This is zero for the mass values (|24p and the redefinition of r mentioned above. This makes it 

possible to choose the conformal gauge u> = h, such that the scalar mode vanishes everywhere. 

The extra conformal gauge symmetry is somewhat reminiscent of the 'partially massless' 

modes that occur in two-derivative Fierz-Pauli theory [la . Ilg] . At the critical value of the 



Higuchi bound Fierz-Pauli theory also develops an extra gauge symmetry 17j j. although not a 
conformal one. The extra gauge symmetry for the higher-derivative theories considered here 
can be thought of as a generalization of the two-derivative partially massless case. 

5 Conserved charges 

We now derive the conserved charges of our theory. They can be calculated via the Abbott- 
Deser method [l8|] , which is an extension of the ADM energy [l9|, |20( to backgrounds with 



constant curvature. In this method the linearized equations of motion E pu are treated as 
an effective energy-momentum tensor. This allows us to compute conserved charges Q^ as 
follows: 

id-l. 



Q"(0 = / d^xV^EFtv. (29) 

JS 

Here £^ is a Killing vector of the background, and £ is a spatial (d — 1) dimensional hyper- 
surface. For instance, the global mass of a solution is then given by Qq for a time-like Killing 
vector. The trick for calculating the conserved charges is to show that the integrand can be 
written as a divergence of a two- form F^ v , 

E^l v = V v F^. (30) 

The integral in (j29[) then reduces to a surface integral at spatial infinity, 



Q"(0= / dS a F^, (31) 

where <9£ is the (d— 2) dimensional boundary of S. For Einstein-Hilbert gravity, whose linear 
equation of motion is simply Q o fo„„ = 0, the two- form is 

f™ = Z p V^K ]p + £ {l y u] h - ^v p h v]p + hp [t y u] i P + ^hVptv 

= ^-°V- ( 32 ) 

Here we have introduced the two-form operator Fe. From the definition above, we have the 
following property. When it acts on a symmetric tensor, it gives a two-form whose divergence 
is the contraction of the Einstein operator with a Killing vector: 

V v F-^oA pu = eQoA pv . (33) 

In our case the linear equation of motion is (|18p . Its general structure is the same as that 
of Einstein gravity, namely a symmetric tensor hit by the Einstein operator. Hence the 
corresponding two-form simply reads 

F liU = ^F i o [j25 + mn oh„ v . (34) 



Like 2jJ,|22j] we restrict to solutions that are asymptotically (A)dS. That is, at spatial infinity 



the vacuum Einstein equations are satisfied: 

0, S 



G pu 



0, R L 



as ^ 



= 0. (35) 

as 



as 

The last equation follows from the fact that the linear Schouten tensor can be decomposed 
as Sn V = Gu V + \~i^9v-yF L . So the terms with Schouten operators in (jM)) are zero in the 
asymptotic region, and all that remains is the product of the squared masses mf. Suppressing 
the dependency on the Killing vector, we obtain 

Q» = ^l[ml (36) 

i=\ 

Thus the conserved charges are equal to those of two-derivative Einstein-Hilbert gravity, up to 
a renormalization factor. In the limit when all extra graviton modes become infinitely heavy 
and decouple, the renormalization factor goes to one by (|10p . Furthermore the conserved 
charges vanish when one of the graviton masses is zero, which is what happens at the critical 
point in four-derivative critical gravity [5L Il2j|. 



6 Graviton energies 

(n) 

In this section we will derive the energies associated to the different graviton modes h^J . 
These modes are annihilated by a single factor of the product in the complete equation of 
motion ()23ap . 

(n-2A-m 2 n )h$ = 0. (37) 

In [23J, y, [5[ the graviton energies were computed by first deriving the Hamiltonian from an 
effective action. The Hamiltonian was then evaluated on-shell for the different graviton modes 
in order to give the energy. The disadvantage of this approach is that one has to make an 
ADM-like split of the indices and variables, and use the Ostrogradsky method to deal with 
the higher derivatives. 

Here we will follow a different route, as outlined in [24], that circumvents these inconve- 
niences. First we compute the energy-momentum tensor by varying the quadratic action (|14[) 
with respect to the background metric: 

V = #• ^ 

The energy is then obtained by integrating this energy-momentum tensor T^ v over a Cauchy 
surface, 

£ = J d^r^nf '■ (39) 

Here n M is the unit normal to £ and £" is a time-like Killing vector. 

For Einstein-Hilbert gravity, the on-shell energy-momentum tensor is 

For deducing the energy-momentum tensor of our theory we need one last identity involving 
the Schouten and Einstein operators. First note from (|19p that when the Schouten and Ein- 
stein operators act on a transverse and traceless tensor, the resulting tensor is also transverse 
and traceless. Furthermore, by equation (|22p . their action on transverse traceless tensors 
gives the same result. This implies that for arbitrary transverse traceless symmetric tensors 
A^y and B^, we have 

B^6 g {S o AfJ) = B^{5- g Q) o A M „ + {SgA,J) Q o B^ . (41) 

Because we will evaluate the energy-momentum tensor on-shell, hu,v is transverse and traceless 
by (|23p . Thus we are allowed to use the above identities in deriving the energy-momentum 
tensor. Lastly, from equations (|22p with (|37p . we have on-shell 

Soh$ =goh$ = -\mlh$. (42) 

Combining the above equations, the energy-momentum tensor becomes 

1 r— 1 r—1 



T (n) = I 
A" 7 r 



\\{mf - m 2 n ) ~ m 2 n J2Yl( m2 j ~ m l) 



=1 3=1 



C"- (43) 



The superscript (n) indicates that it is evaluated on-shell for the mode njjj . When m n = 0, 
only the first product contributes to the energy-momentum tensor. For non-zero masses rrii, 
the first product vanishes and only one term in the sum is non-zero. This gives the following 
energies for the massless and massive gravitons: 

r— 1 
f (0) =+ ^HjJ m 2 (44a) 

3=1 

r— 1 
g(i) = _^l m | JJ( m J _ m 2). (44b) 

T i=i 

The energy of the massless graviton h^J has the same overall factor as the conserved charge 

/.a 

Q^ ()36p . So it seems that the massive gravitons h)iv do not contribute to the conserved charge. 

This is to be expected, as the massive gravitons fall of too fast towards spatial infinity to 

contribute to the surface integral (|31|) in the asymptotic region. 

For four-derivative theories (r = 2), the above energies reduce to £' ' = —S^ 1 ', which 
matches with the energies found in [J, |5(. So for rank two the only way to obtain energies 
with the same sign is to set the mass mi to zero. Both energies are then zero, and the 
conserved charge also vanishes, rendering the theory trivial. 

In the six-derivative case, r = 3, this zero-energy problem does not occur. The graviton 
energies (|44p then namely read 

5(0) = ^ELmlml (45a) 

r 

< r(D = ^i m 2 (m 2_ m 2 )) (45b) 

r 

£ m = ^i m 2 {m 2_ m 2 ) _ (45c) 

r 



These energies are plotted in Figure 1 There are two distinct points where they have same 
sign: either when mf = m 2 , or when m\ = V m 2 , = 0. In the last critical point the conserved 
is zero, whereas it can be positive in the first. 

For yet higher rank theories there are even more critical points. However, it will never 
be possible to have the same sign for all energies without degeneracies in the masses. If we 
namely arrange the masses by size, 

m\< ... <m 2 i < ... < m 2 _ l5 (46) 

the sign of the energies (|44p alternates: 

sgn (s®) = - sgn (£ (i+1 A . (47) 



Thus the situation with all masses different (j4"6j) leads to ghosts; to avoid this we need at least 
some degeneracy of the masses. 

But whenever there is a mass m n with multiplicity \i n > 1, a so-called log-mode nf^f 



appears j8|]. These log- modes are annihilated not by a single, but by multiple factors of the 
product in the equations of motion, 

(a-2A-ml) p h^ = 0. (48) 
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Figure 1: Graviton energies for r = 3 in polar coordinates. Here 9 = tan l ( ^f- J is the angle in the 

(m^mf) plane. The masses are given by m\ = 2cos# and m| = 2sin6*. There arc three 
points where all energies are non- negative: m\ = mf, m\ = 0, or m| = 0. 



The label p can take the values p = 2, . . . ,p n , as p = 1 simply gives the non-logarithmic 

(n) 

raviton mode h^J . From the four-derivative case the log- modes are expected to be ghosts 
and, if possible, need to be truncated out in order to restore unitarity. 



7 Conclusions 



In this paper we have studied gravities of arbitrary rank, meaning they propagate any number 
of gravitons on (A)dS backgrounds. Besides from giving a quadratic and a non-linear action, 
we have calculated the conserved charges and the graviton energies. From the energies we 
deduce that there will be ghosts unless the masses have critical values. At these critical points 
some of the gravitons have degenerate mass. But as mass degeneracies lead to logarithmic 
graviton modes, the untruncated theory will never be unitary. By truncating the log-modes 
by imposing appropriate boundary conditions one could obtain a unitary sub-sector of the 
theory. We leave the exact form of both the higher-rank log-modes and boundary conditions 
to future study. 

When the rank is two, there is only one critical point, and all the energies vanish [4.l5l.ll2|. 
One can interpret the triviality of this theoryas being to due to the proposed equivalence of 
Einstein gravity and conformal gravity 25|, |26[ . We have shown that for higher rank theories 
there are critical points where the conserved charges and graviton energies do not vanish. But 
in the fully degenerate case where all the graviton masses are zero, the theory will always be 
empty. Like the proposed equivalence of Einstein and conformal gravity, this 'emptiness' of 
higher rank theories could in principle be used to construct a chain of equivalence relations 
between gravity theories of different rank. 
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A Conventions 

We use the 'mostly plus' metric signature (— ,+, . . . .+). The conventions for the Riemann 
tensor are the default of the xAct software package [22], which in turn follows Wald's conven- 
tions hi 



[v pi » U}-*- p — -ttpvp 1<ji 



R 



i m/ 



R p 

1 l ppv ■ 



(49) 



Barred objects are background quantities (i.e. g denotes the background metric). AdS and 
dS backgrounds are chosen as follows: 



Hfivpcr — A \9pp9va 9pa9vp) 



R 



pv 



(d - 1)AS 



pvi 



R = d{d-1)A. 
Perturbations around these backgrounds are defined as 

9pv = 9pv ~r 9pu = 9pv ~r tlpv 



(50a) 
(50b) 
(50c) 



(51) 



The superscript L indicates linear perturbations. Thus the linear perturbation of the metric 
is given by h^ u . 

The cosmological Einstein tensor is the usual Einstein tensor plus a term proportional to 
the cosmological constant, such that it vanishes on the above backgrounds: 



^ pv 
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1 



pv 
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pv 



The cosmological Schouten tensor is defined similarly, 



S: 



S: 



pu 

A 
pv 



R 



1 



pv 



R 



-flpv \R-[d- 2) (d - 1) Aj , (52a) 

0. (52b) 

-<!,,„ j— { +(d-2)A , (53a) 

0. (53b) 

The Schouten tensor is usually given with an additional overall factor -g—n. However, for 
our purposes the above definition is more convenient. In the main text the superscripts A 
are dropped from the cosmological Einstein and Schouten tensors. Thus by G^ and S^ v we 
always mean their cosmological versions. 

For completeness, we give the linear perturbations of the cosmological Einstein and 
Schouten tensors: 
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pi* 



S, 



pv 



1 (d - l)A/v 

(d-l)AV 



pv K" 1 *-j±**«fu> 
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pv 



-j9pv R 
1 

2{d-l 



-9pvR , 
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pv 



Vp^ipKf 
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-ah 



l. 



pv 



-V„V„/i 
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R L 



VpVah pa -Uh-{d- l)Ah. 



(54a) 
(54b) 

(55a) 
(55b) 
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